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We propose a measure of non-Gaussianity for quantum states of a system of n oscillator modes. 
Our measure is based on the quasi-probability Q{a), a £ C". Since any measure of non-Gaussianity 
is necessarily an attempt at making a quantitative statement on the departure of the shape of the 
Q function from Gaussian, any good measure of non-Gaussianity should be invariant under trans- 
formations which do not alter the shape of the Q functions, namely displacements, passage through 
passive linear systems, and uniform scaling of all the phase space variables: Q{a) — > \^"Q{Xa). 
Our measure which meets this 'shape criterion' is computed for a few families of states, and the 
results are contrasted with existing measures of non-Gaussianity. The shape criterion implies, in 
particular, that the non-Gaussianity of the photon-added thermal states should be independent of 
temperature. 
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PACS numbers: 03.67.-a, 03.65.-w, 02.50.Cw 



I. INTRODUCTION 

Quantum information theory of continuous variable 
systems has been actively pursued in recent years, es- 
pecially in the context of Gaussian states [U, 0, [1] ■ Such 
states are the ones which occur naturally in most experi- 
mental situations, particularly in quantum optics. While 
these states live in an infinite-dimensional Hilbert space, 
they are remarkably easy to handle since they are fully 
described by their covariance matrix (and first moments). 
Further, their evolution under quadratic Hamiltonians 
is easily cast in the language of symplectic groups and 
(classical) phase space [j, [M 13 ■ The fundamental pro- 
tocol of quantum teleportation has been achieved using 
these states 0, II]- However, there are situations wherein 
one deals with (nonclassical) non-Gaussian resources to 
generate entanglement [1, [M [H [H, [H [11 d, [Hi- 
They arise naturally in nonlinear evolutions like passage 
through a Kerr medium [13, [3 • 

It has been shown recently that teleportation fideli- 
ties can be improved with the use of non-Gaussian re- 
sources [l^ . It is thus important that one is able to 
quantify the non-Gaussianity of such resources. Effort 
in this direction has been initiated in some recent publi- 
cations [13, HI . 

From the perspective of classical probability theory, 
Gaussian distributions are those probability distributions 
which are completely specified by their first and sec- 
ond moments; all their higher-order moments are de- 
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termined by these lower-order moments. Non-Gaussian 
probabilities do not enjoy this special property. An eas- 
ier, and possibly more effective, way to distinguish the 
two is through cumulants: every non-vanishing cumulant 
of order greater than two serves as an indicator of non- 
Gaussianity of the probability distribution under consid- 
eration [H, . 

The purpose of any good measure of non-Gaussianity 
in the context of classical probability theory is thus to 
capture the essence of the non- vanishing higher-order cu- 
mulants. A non-Gaussianity measure should thus mani- 
festly depend on the higher-order cumulants. Yet another 
desirable feature one would like to have is invariance of 
the measure under scaling. Ultimately, non-Gaussianity 
measure is a quantitative statement of the departure of 
the shape of a probability distribution from Gaussian. But 
uniform scaling of all the variables of a probability distri- 
bution does not alter the 'shape ' of the distribution, and 
hence it should not affect its non-Gaussianity. 

The notion of non-Gaussianity can be extended to a 
quantum mechanical state through its definition on the 
associated Q function, a member of the one-parameter 
family of s-ordered quasi-probabilities [l^l . That this 
is an appropriate route is endorsed by the fact that 
the Marcinkiewicz theorem [see below] holds for the s- 
ordered quasi-probabilities as well. It turns out that 
the cumulants of order greater than 2 for the various s- 
ordered quasi-probabilities corresponding to a fixed state 
p are independent of s, indicating that the higher or- 
der cumulants are intrinsic to the state. Moreover, all 
higher-order cumulants of order greater than 2 vanish 
identically for Gaussian states. Thus any non- vanishing 
higher-order cumulant of the quasi-probability indicates 
non-Gaussianity of the state, and this conclusion is inde- 
pendent of the ordering parameter s. 
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The above considerations will suggest that any good 
measure of non-Gaussianity relevant in the context of 
classical probability theory can, with suitable modifi- 
cation, lead to a good measure of non-Gaussianity of 
quantum mechanical states, provided a state is identi- 
fied through its Q function ( For a brief review on such 
measures in classical probability theory, see [23| ). The 
purpose of such a quantum measure would be to capture 
the essence of the non- vanishing higher-order cumulants 
of the Q function associated with the state. And invari- 
ance of the measure under an overall scaling of the Q 
function is a desirable feature worth insisting on. 

The purpose of this paper is to motivate and present 
such a measure of non-Gaussianity of quantum states. 
Our measure is based on the Wehrl entropy [25*1 , the 
quantum analogue of differential entropy |26i] well-known 
from the context of classical information theory of con- 
tinuous variables [Differential entropy itself is a general- 
isation of Shannon entropy from discrete to continuous 
variables] . 

The photon-added thermal states [131 play a key role in 
our considerations. These nonclassical states have been 
generated experimentally [H, [2^, [s^, [sH . Their special 
importance to the present work arises from the fact that 
the Q functions of these states are scaled versions of those 
of the Fock states, and therefore one will expect any good 
measure of non-Gaussianity to return the same values for 
both classes of states. 

The plan of the paper is as follows. In Section II we 
briefly introduce the definition of moments and cumu- 
lants, and recall two well-known theorems in the context 
of these notions. The one-parameter family of s-ordered 
quasi-probabilities corresponding to quantum density op- 
erators is briefly considered in Section III, with particu- 
lar emphasis on the Q function, and in Section IV we 
review the relationship between differential entropy and 
the KuUback-Leibler distance of classical probability the- 
ory. As a final item of preparation, we review briefly in 
Section IV the Wehrl entropy [1^ and some of its prop- 
erties. With these preparations, we introduce in Section 

VI our non-Gaussianity measure and explore some of its 
more important properties, including its invariance under 
uniform scaling of the underlying phase space. In Section 

VII we evaluate this measure for three families of quan- 
tum states, and in Section VIII we compare our measure 
with two other measures of non-Gaussianity available in 
the literature. The paper concludes in Section IX with 
some additional remarks. 



II. MOMENTS AND CUMULANTS 

For a multivariate probability distribution 7^ (a;), where 
X = {xi,X2, ■ ■ ■ ,Xn) € TZ^ , the characteristic function 
x(0: G "'^"i is given by the Fourier transform of 
V{x) O: 

X(^) = / d''xV{x)eyi-p[i£,-x] 



<.r.r-<-)=/^"..r.r-<-n.).(2.i) 

It follows from the invertibility of Fourier transformation 
that the characteristic function retains all the informa- 
tion contained in the probability distribution. The char- 
acteristic function is often called the moment generating 
function, since one obtains from it all the moments of the 
underlying probability distribution through this compact 
expression: 

Another equivalent description of a probability dis- 
tribution is through the cumulant generating function. 
This is defined through the logarithm of the characteris- 
tic function 

= E (n^)7™......,.„,(2.3) 

iTiim2---m„ \fe=l ' / 

or, equivalently, through 

xiO = exp(r(e)) . (2.4) 

From Eq. p.3p . it is easy to see that the cumulants 
7mi,m2, ■■■,'"„ can be expressed as 

7mi,m2,--,mfe = ( d{i£,k)""' ) ^^^^^ ' ^^'^^ 

Thus, the cumulants are related to r(-) in precisely the 
same way as the moments are related to x(-)- The set 
of all moments {x'^^x'!^^ ■ ■ ■ x™") gives a complete char- 
acterisation of a probability distribution V{x), and the 
same is true of the set of all cumulants "fmi,m2.---,m„ as 
well. Indeed, one can describe one set in terms of the 
other [22, 32, 33]. 

With these notations and definitions on hand, we now 
recall two important results from classical probability 
theory. 

Theorem 1 The cumulant generating function of a 
Gaussian probability distribution in n variables is a 
multinomial of degree equal to 2 fjl /. 

Theorem 2 (Marcinkiewicz Theorem). If the cumulant 
generating function of a (normalised) function in n vari- 
ables is a multinomial of finite degree greater than 2, 
then the function will not be point wise non-negative, and 
hence will fail to be a probability distribution ]3A . \ 3aj . 
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Theorem 1 is a statement of the fact that a Gaus- 
sian probabihty is fully determined by its moments of 
order < 2 ; all the higher-order cumulants are identically 
zero for a Gaussian probability. Theorem 2 is a much 
stronger statement. It implies that any true probability 
distribution other than the Gaussian distribution has a 
cumulant generating function which cannot truncate at 
any (finite) order. That is, a non-Gaussian probability 
distribution has non-vanishing cumulants of arbitrarily 
high order. We note in passing that non-vanishing cu- 
mulants of order greater than 2 serve as indicators of the 
non-Gaussianity of the underlying probability. 

III. QUASI-PROBABILITIES AND THE Q 
FUNCTION 

A state of a quantum mechanical system specified by 
density operator p can be faithfully described by any 
member of the one-parameter family of s-ordered quasi- 
probability distributions — 1 < s < 1 In other words, 
an s-ordered quasi-probability captures all the informa- 
tion present in the density operator p. However, it is not 
a genuine probability distribution in general; in particu- 
lar, it is not point wise non-negative. The prefix quasi 
underscores precisely this aspect. Nevertheless, the s- 
ordered family of quasi-probability distributions gives us 
a framework wherein one could give a phase space de- 
scription of quantum mechanical systems in the language 
of classical probability theory. 

For a quantum state describing the radiation field of n 
modes {n oscillators) the characteristic function of the s- 
ordered quasi-probability, for any — 1 < s < 1, is defined 
through [2J] 

Xp(C;s)=exp[^|enTr(pi^(0), (3.1) 

where ^ = (Ci; Ca, • ' • j ?«) G C", and D{^) is the n-mode 
(phase space) displacement operator: 

D{0 = exp[^(e,a]-r«,)]. (3.2) 

3 

The s-ordered quasi-probability itself is just the Fourier 
transform of this characteristic function XpiCi s)' 

Wpia; s) = /'exp[5](a*C.- - a,^*) ]xpiC, s) J] d^Cj-^ 

•' 3 3 

Here and a}- are the annihilation and creation opera- 
tors of the jth mode, aj represents the (c-number) phase 
space variables qj^Pj corresponding to the jth mode 
through — {qj -\- ip.j)/\/2, and a = (ai, q;2, • • • , a„) G 
C". The particular cases s = —1,0,1 correspond, re- 
spectively, to the better known Q function, the Wigner 
function, and the P function. 

The Q function corresponding to a density operator p 
has a particularly simple expression in terms of coherent 



state projections: 

Q{a) = {a\p\a), a € C" . (3.4) 

It may be noted that the Q function is manifestly non- 
negative for all a e C". 

Reality of Wp{a;s) is equivalent to hermiticity of the 
density operator p, and the fact that p is of unit trace 
faithfully transcribes to 

^ J Wpia;s)d^a^l. (3.5) 

While these two properties hold for every s-ordered quasi- 
probability, point wise non negativity for all states is a 
distinction which applies to the Q function alone. In 
other words, the Q function is a genuine probability dis- 
tribution; every other Wp{a] s) is only a quasi-probability. 
Gaussian pure states are the only pure states for which 
the Wigner function is a classical probability [IBl ; in the 
case of P function, the coherent states are the only pure 
state with this property. 

However, not every probability distribution is a Q func- 
tion. This is evident, for instance, from the obvious fact 
that Q{a) < 1, Va G C". 

The next result captures, in a concise form, the man- 
ner in which members of the one-parameter family of 
s-ordered quasi-probabilities Wp{a; s) differ from one an- 
other for a given state p. 

Theorem 3 Only the second order cumulants of the 
quasi-probability of a given state depend on the order pa- 
rameter s; all the other cumulants are independent of the 
quasi-probability under consideration. 

This result is already familiar in the case of a single- 
mode radiation field [17|. But the proof is, as outlined 
below, immediate in the multi-mode case as well. The 
characteristic functions of a state p for two different val- 
ues of the order parameter si and S2 are obviously related 
in the following manner 24]: 

Xp(C; si) = exp (^{si - S2)|^|^) Xp(C; S2) ■ (3.6) 

On taking logarithm of both sides to obtain the corre- 
sponding cumulant generating functions we have 

logXp(C;si) = (si -s2)|^|VlogXp(?;s2). 

That is, 

rp(C;si) = (si-s2)|e|Vrp(C;s2). (3.7) 

Thus the cumulant generating function for different s- 
ordered quasi-probabilities differ only in second order, 
completing proof of the theorem. 

In these equations stands, as usual, for l^jl^- 
As an immediate consequence of this theorem we have 

Theorem 4 For no quasi-probability can the cumulant 
generating function be a multinomial of finite order > 2. 
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Proof: Since the Q function, for every state p, is 
a genuine probability distribution, it follows from the 
Marcinkiewicz theorem that the cumulant generating 
function of Q cannot be a multinomial of finite order 
> 2. Since the different s-ordcrcd quasi-probabilities dif- 
fer only in second-order cumulants, this conclusion holds 
for all s-ordered quasi-probabilities, thus proving the the- 
orem. 

We conclude this Section with the following remarks. 
The above considerations show that quasi-probabilities 
fail to be true probabilities only in this limited sense: 
they differ from genuine probabilities only in cumulants 
of order two. The distributions, however, can be quite 
different from classical probabilities, particularly for s > 
0, and they can become as subtle as Fourier transform of 
exp [cry^], a > 0, a Gaussian with the wrong signature 
for the variance. 

Since the higher-order cumulants, which should play 
an essential role in any reasonable definition of non- 
Gaussianity measure, do not depend on the value of the 
parameter s, they may be viewed as attributes intrinsic to 
the state under consideration; we may therefore use any 
convenient quasi-probability to capture their essence. 

IV. DIFFERENTIAL ENTROPY AND THE 
KULLBACK-LEIBLER DISTANCE 

The role of Shannon entropy of probability distribu- 
tions over discrete random variables is taken over by 
differential entropy in the case of continuous variables. 
Given a multivariate probability distribution V{x) in n 
variables {xi, X2, ■ ■ ■ , Xn) & 7?,", the associated differen- 
tial entropy H{V{x)) is defined by [l^ 

H{V{x)) = ~ J d''xV{x)\ogV{x) . (4.1) 

But unlike the Shannon entropy, the differential entropy 
can be negative. This is manifest, for instance, for uni- 
form distribution over a region of less than unit volume 
in 7^". 

Among all the probability distributions with a fixed 
set of first and second moments, the Gaussian probability 
distribution has the maximum differential entropy ^26|. 
This fact may be used to modify differential entropy to 
result in a non-negative quantity 

J(rix)) = HiVG{x))-HiV{x)). (4.2) 

Here Vg {x) is the Gaussian probability distribution with 
the same first and second moments as the given proba- 
bility distribution V{x). 

It may be recalled that KuUback-Leibler distance be- 
tween two probabilities Pi (x) and P2 (x) is defined as the 
difference of their differential entropies [2^ : 

Siriix)\\r2ix)) = H{V2{x)) - H{Vi{x)) 

^ - J Vi{x)log{Vi{x))d"x 



+ J V2{x)\og{V2{x))d"x. (4.3) 

Thus J{V{x)) can be regarded as the KuUback-Leibler 
distance between the given probability P{x) and the as- 
sociated Gaussian distribution Pq{x): 

J{nx))=S{VG{x)\\V{x)). (4.4) 

J{V{x)) is sometimes known by the name negentropy. 

V. WEHRL ENTROPY 

Wehrl entropy [2^ [s^l niay be viewed as the extension 
of differential entropy to the quantum mechanical con- 
text, but the Wehrl entropy has interesting properties 
which distinguish it from differential entropy. The dis- 
tinction arises from the fact that while every Q function 
certainly qualifies to be a classical probability distribu- 
tion, every classical probability is not a Q function. The 
uncertainty principle has a fundamental role to play in 
this aspect [25| . The potential use of Wehrl entropy as a 
measure of the 'coherent' component of a state has been 
discussed in Ref [s^ ■ And its possible role in defining an 
entanglement measure has also been explored [s^, |40| . 

For a state p describing n modes of radiation field, the 
Wehrl entropy is defined as 

Hw{p) = / \{d^a,Qp{a)\ogQp{a), (5.1) 

where Qp{a) is the Q function corresponding to p. This 
definition may be compared with that of differential en- 
tropy; the role of 'P{x) in differential entropy is played 
by Qp{a) in Wehrl entropy. 

However, in contradistinction to differential entropy, 
the Wehrl entropy is always positive. This is an im- 
mediate consequence of the fact that the Q function is 
bounded from above by unity. It turns out that the Wehrl 
entropy is always greater than or equal to unity [41] ; in- 
deed, it attains its least value of unity for the coherent 
states and only for these states. This property can be 
thought of as a manifestation of the uncertainty prin- 
ciple, which the coherent states saturate. Further, the 
Wehrl entropy is always greater than the von Neumann 
entropy [1^1 : 

Hw{p)>S{p)^-TY{p\ogp). (5.2) 

While the von Neumann entropy is zero for pure states, 
we have just noted that the Wehrl entropy Hw{p) is 
greater than or equal to unity for all states. Several 
aspects of the Wehrl entropy have been explored in 
Ref. [H. 

VI. A NON-GAUSSIANITY MEASURE FOR 
QUANTUM STATES 

As is well-known, a quantum state p is said to be Gaus- 
sian iff the associated Wigner distribution is Gaussian. 
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This will suggest that the non-Gaussianity of a state is 
coded into the non-vanishing cumulants of order > 2 of 
the Wigner function. Since the Wigner and Q functions 
are related by convolution by a Gaussian, the Q func- 
tion of a state is Gaussian iff the Wigner function is, and 
the non-Gaussianity should thus be found coded in the 
higher-order cumulants of the Q function as well. The 
consistency of these statements is ensured by the fact that 
the higher-order cumulants are the same for the Wigner 
and the Q functions [Indeed, as we have shown earlier, 
the higher-order cumulants are intrinsic to the state, and 
hence are the same for all s-ordered quasi- probabilities] . 

Non-Gaussianity can thus be described using either the 
Wigner function or the Q function. The fact that the Q 
function is everywhere non-negative, rendering it a gen- 
uine probability in the classical sense, makes it our pre- 
ferred choice. We employ therefore the Wehrl entropy to 
capture the essence of the higher-order cumulants. 

Given a state p, our measure of non-Gaussianity M{p) 
is defined as the difference of two Wehrl entropies: 

Af{p) = Hw{pg)-Hw{P). (6.1) 

Here Hw (p) is the Wehrl entropy of the given state p and 
Hw{pg) is the Wehrl entropy of the Gaussian state pa 
that has the same first and second moments as p. Since 
M{p) measures the departure of the Wehrl entropy of p 
from that of its Gaussian partner pG^ it can be viewed as 
a quantum KuUback-Leibler distance. N{p) could also 
be viewed as a relative Wehrl entropy. But we prefer to 
call it simply a non-Gaussianity measure. 

This measure of non-Gaussianity enjoys several inter- 
esting properties. We will now list some of them: 

(i) N{p) > 0, equality holding iff p is Gaussian. 
Proof: This is a restatement of the fact that the Wehrl 
entropy of a Gaussian state is greater than that of all 
states with the same first and second moments as the 
Gaussian. 

(ii) A/'(/5) is invariant under phase space displacements: 

Nip)=N{D{OpD0). (6.2) 

Proof: Let D{£^) p D{£^y be denoted, for brevity, by p'. 
The Q function of p ' is related to that of p in this simple 
manner: 

Qp'{a) = Qp{a-0. (6.3) 

That is, displacement i?(^) acts as a rigid translation in 
phase space '2^ |4ll, \^ . Thus it has no effect on the 
Wehrl entropy of any state, and hence leaves Af{p) in- 
variant for every state. 

(iii) M{p) is invariant under passage through any pas- 
sive linear system. 

Proof: A passive linear system is represented by a n x n 
unitary matrix U . It maps a coherent state |a) into a new 
coherent state \a') = | f/ a) [H |4l|. Ell , where a G C" is 
to be viewed as a column vector. Let lAu be the unitary 
operator in the n-mode Hilbert space which represents 



the passive linear system labelled by the matrix U . Let us 
denote by p' the transformed state Uu pU^ at the output 
of this passive system. Then the output Q function is 
related to the input Q function in this manner: 

Qp, [a) = Qp{U-^a) = Qp{U^a) . (6.4) 

That is, the action of a passive linear system is a rigid 
S0{2n) rotation in the 27i-dimensional phase space. It 
follows immediately that this transformation does not 
change the Wehrl entropy of any state, and hence does 
not affect M{p). 

Remark: While in the single-mode case of two- 
dimensional phase space all proper rotations are canon- 
ical transformations, this is not true in the multi-mode 
case. That is, Sp(2n, 7^) n S0(2n) is a proper subgroup 
of S0(2n) isomorphic to U(n), the n^-parameter group of 
n X n unitary matrices, whereas S0(2n) is a much larger 
(2n^ — n)-parameter group [i^ • Only those phase space 
rotations which are elements of this intersection act as 
unitary transformations in the Hilbert space of n oscilla- 
tors. 

(iv) M{p) is invariant under a uniform phase space scal- 
ing A defined at the level of the Q function in the following 
manner: 

A : Q{a) -> Q'{a) = A2"Q(Aa) . (6.5) 

Proof: Under this uniform phase space scaling of the Q 
function, the Wehrl entropy changes by a simple additive 
part that is independent of the state: 

If " 
Hw{p) = -— / Q{a)logQ{a)Y[d^aj 

^ -— / A2"Q(Aa)log(A2"Q(Aa))n'^'«, 

= Hw{p)-2n\og\. (6.6) 

Note that in arriving at the last equation we have made 
a change of variables in the integral and made use of 
the normalisation of the Q function. Now it trivially 
follows from this result that M{p), being a difference of 
two Wehrl entropies, remains invariant. 
Remark: While the above conclusion holds mathemati- 
cally for all A > 0, the scaled Q function fails to be a phys- 
ical Q function if A > 1. Therefore we restrict this scale 
parameter to the physically relevant range < A < 1. 

(v) M{p) is additive on tensor product states: 

^f{Pl®p2) = ^f{pl)+^{p2)■ (6.7) 

Proof: Under tensor product the Q functions go as prod- 
uct probabilities by definition. This is true of their asso- 
ciated Gaussian probabilities as well. 
Corollary : For a bipartite state of the form p = Pa® PG^ 
where pG is a Gaussian state 

^^{p)^^^{pa®PG)=^f{Pa)■ (6-8) 
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Proof: From (v) we have 

J\f{p) = U{pa ® Pg) = J^{Pa) + J^{Pg) ■ 

and from (i) 

U{pa)+^f{pG)=M{pa). (6.9) 



Proposition: For a bipartite state of the form pout = 
i^u {pa 'S) \a){a\)U^, where U represents a passive hnear 
system and la) is a coherent state, we have 

^fiPout)^^fiPa)■ (6.10) 

Proof: From (iii) we have 

AA(/5out) = J^iUu {Pa ® \a){a\)Ulj)^JVipa ^\a){a\). 

We have from (v) 

Af{pa(g>\a){a\) =Afipa)+Afi\a){a\). 

Since the coherent state |a) is Gaussian, we have from 
(i) 

Af(pa)+Af{\a){a\)=Af{pa). (6.11) 

This result is useful in evaluating the non-Gaussianity 
of bipartite states produced by the action of beam split- 
ters, as we shall illustrate in the next Section. 



A. Shape criterion for good measure of 
non-Gaussianity 

Properties (ii), (iii), and (iv) deal with transforma- 
tions which do not change the shape of the Q functions. 
Since non-Gaussianity is a quantitative statement regard- 
ing the departure of the shape of the Q function from 
Gaussian, it will appear that any good measure of non- 
Gaussianity should return the same value for all states 
connected by these transformations. In particular, two 
quantum states whose Q functions are related by a uni- 
form scaling of all the phase space coordinates should be 
assigned the same amount of non-Gaussianity. We will 
call this the shape criterion, and we have seen that our 
measure Af{p) meets this requirement. 



VII. EXAMPLES 

In this Section we evaluate our non-Gaussianity mea- 
sure Af{p) for three families of states, namely the Fock 
states, the photon-added thermal states, and the phase- 
averaged coherent states of a single-mode of radiation. 
While the first two families consist of nonclassical states, 
the third one is a family of classical states. 



A. Photon number states 

The Q function of the Fock state (energy eigenstate) 
p = \m){m\ of the oscillator is given by the phase space 
distribution 



?|m>(a) = -^^-exp(-|a|^) , 



(7.1) 



whose only non-vanishing moment of order < 2 is 
(|a|^) — Tr(paat) = m + 1. The phase space aver- 
age (|a| ) is with respect to the probability distribu- 
tion Q\m)ict) and, by definition, it equals the (quan- 
tum) expectation value of the associated anti-normally 
ordered operator aa^. The Gaussian state which has 
the same moments of order < 2 as p\m) = \m){m\ 
is clearly the thermal state with mean photon number 
m. The Q function of such a thermal 



state Pg is given by 



QGia) 



1 



exp 



1 



(7.2) 



The Wehrl entropy corresponding to pG is easily com- 
puted: 



HwiPG) = l + log(l + (n)) 
= l + log(l + m). 



(7.3) 



The Wehrl entropy of the photon number state p 
\m) (ml is 



Hw{p\m)) 



c^^a'9|m)(a)log(9|,„)(a) . (7.4) 



This can be computed explicitly by going to the polar 
coordinates, and one obtains [38] 



Hw{p\m)) = I + m + \ogm\ — mip{m + 1), 



lib -j^ 

■0(^ + 1) = X! ^ ^ 
fc=i 



(7.5) 



where 'ip{m) is the digamma function, and 7 = 0.5772 . • . 
is the Euler constant. Hence the non-Gaussianity of the 
photon number state p ~ \m){m\ is 

^{P\m)) = Hw{Pg) - Hw{p\m)), 

— ln(m + 1) — m — logm! + mil>{m -\- 1).(7.6) 

In Fig.lJT]) we have plotted this non-Gaussianity as a 
function of the photon number m. It is clear that the 
non-Gaussianity of \m) increases monotonically with the 
photon number m, and goes to 00 as to tends to 00. That 
this was to be expected can be seen as follows. For large 
TO values '0(to-|-1) ^ ln(TO+l), and \ogm\ ~ TOlogm — m, 
and hence N{p\m)) ~ log(TO -I- 1). We shall be returning 
to this result in the next Section. 

Now consider a bipartite state of two modes with one 
mode in the Fock state and the other in the vacuum. 
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Non-Gaussianity of this product state is the same as that 
of the Fock state, and this follows from Eq. (|6.8p . Let 
this bipartite state be passed through a beam splitter. 
The state at the output will be entangled due to the 
nonclassicality of the Fock state 0, , but in view of 
Eq. ()6.10p . this two- mode state will have the same non- 
Gaussianity as the original single-mode Fock state. 



B. Photon-added thermal states 

In this subsection we evaluate the non-Gaussianity of 
the photon-added thermal state (PATS) 0\. The PATS 
is defined through 



(7.7) 



where C is the normalisation constant which ensures 
Tr (p) — 1, and pth is the thermal state given by 



Pth 



(1 — x) ^ a;'^|fc)(fc| ; x = exp 



Tl = 



kf 



(7- 



One can alternatively define the PATS through paramet- 
ric differentiation: 



m+l ,„i 



k=Q 



X^\k){k\ 



(7.9) 



PATS are thus parametrised by two parameters: < 2: < 
1, and TO = 0, 1, 2, • • • . The limit a; — > corresponds to 
Fock states, and the limit m — > corresponds to thermal 
states. 

We may note that PATS (with to > 1) is nonclassical 
for all values of a; Q . Indeed, it violates a three-term 
classicality condition . 

The Q function of PATS can be easily calculated and 
is given by 

Qi^Tk^) = ^^^^] |a|""exp[-(l - x)\af].{7.10) 

It is evident that the Q function of the PATS is a scaled 
version of the Q function of the Fock state: 

QpATl(«)=^'Q|m)(A«): A = VT^. (7.11) 

Since our measure of non-Gaussianity respects the shape 
criterion put forward in the previous Section, it is imme- 
diate that the non-Gaussianity of the PATS is the same 
as that of the photon number state: 



ln(TO + 1) — TO — logTO! + m'ijj{m + 1) 
^(P|™)). (7.12) 



It is worth emphasising here that the PATS is a spe- 
cial state with regard to the question of verifying whether 
a given measure of non-Gaussianity is a good measure, 
i.e., whether it satisfies the shape criterion. The test is as 




FIG. 1: Variation of A/'(p) with number of photons m for the 
Fock state p = \m){m\. 



a. 




FIG. 2: Variation of 7V(p) with energy |/3^| for the phase- 
averaged coherent state. 



simple as checking whether the measure in question eval- 
uated for the PATS is independent of the temperature 
parameter x or not. 

Finally we consider, as in the previous Subsection, a 
bipartite state of two modes, with one mode in the PATS 
PpATS ^^'^ other in the vacuum state. Let us pass this 
two-mode state through a beam splitter. That the state 
at the output of the beam splitter is entangled follows 
from the nonclassicality of the PATS [9|, ilOj . It follows 
from Eq. (|6.10p that non-Gaussianity of this entangled 
state is the same as that of the PATS, and hence is fully 
determined by m. 

We have already noted that PATS violates a three- 
term classicality condition. This implies that the output 
state is entangled and, moreover, that it is single-copy 
distillable[i|. 



C. Phase-averaged coherent states 

As our final example, we evaluate the non-Gaussianity 
for the phase-averaged coherent states. Given a coherent 
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FIG. 3: Variation of Si{p) as a function of the Boltzmann 
parameter x for the photon-added thermal state. 



FIG. 4: Variation of 52 (p) as a function of the Boltzmann 
parameter x for the photon-added thermal state. 



state 1/3) its phase- averaged version is 



VIII. COMPARISON WITH OTHER 
MEASURES 



P\P\ 



de 

2^ 



eyiY)[—iOa)a] |/3)(/3| exp[i6'a^a] 



exp(-|/?|2)5]^|n)H 
^ — ' ni 



(7.13) 



n=0 



Since /5|^| is a convex sum of Fock states, its Q function 
is a corresponding convex sum: 



Ql'^l(a) = exp[-(|a|- 
= cxp[-(|q;|- 



,2,, \a\ 



n=0 

^)]/o(2|a||/3|), (7.14) 



where /o(-) is the modified Bessel function of integral or- 
der zero. The only non-zero moment of order < 2 is 
(|z|^) = TT{p\fs\aa^) = 1 + |/3|^. The associated Gaussian 

probability QQ^{a) that has the same first and second 
moments is thus the thermal state with average pho- 
ton number (n) = As we have shown earlier in 
Eq. (|7.3p . the Wehrl entropy of this Gaussian state is 

Hw{Pg) = 1 + log(l + 1^1^)- To compute the Wehrl 
entropy corresponding to the original phase-averaged co- 
herent state, however, we resort to numerical evaluation. 
In Fig. ([2|) we present the non-Gaussianity of /5|^| as a 
function of the energy of the state. It is seen to be 
a monotone increasing function of 

Note that the phase-averaged coherent states are clas- 
sical since they are, by definition, convex sums of coher- 
ent states. Thus if a bipartite state consisting of a phase- 
averaged coherent state in one mode and vacuum in the 
other is passed through a beam splitter, the two-mode 
mixed state at the output will remain separable (since 
the phase-averaged coherent state is classical ^ ) , with 
the same non-Gaussianity as the original phase-averaged 
coherent state. 



In this Section we compare our non-Gaussianity mea- 
sure M{p) with two non-Gaussianity measures which 
have been proposed recently. 



A. Measure based on Hilbert-Schmidt distance 

Genoni et al [20| , have proposed a non-Gaussianity 
measure based on the Hilbert-Schmidt distance. They 
define non-Gaussianity of a state p as 



Slip) 



2Tr(/52) 



(8.1) 



where f is the Gaussian state with the same first and 
second moments as p. Let us compare this measure with 
ours in the specific case of the PATS /5p™^s- ^^S- ® 
we plot Si (/5p™^s ) ^ function of the Boltzmann param- 
eter X, for fixed value of m = 1. It is seen that ^i(Ppats)' 
for m = 1, is not a constant but varies with the tem- 
perature parameter x. This shows that this measure of 
Genoni et al does not satisfy our shape criterion. 

Another interesting difference appears when one com- 
pares our measure A/'(/5) with Si{p) in the case of the pho- 
ton number states p = \m){m\. As we have shown ear- 
lier [see Fig. ([1])], our measure monotonically increases 
with the photon number m and tends to infinity as m 
tends to infinity. In contrast, as Genoni et al have shown 
and emphasised [2^1 , their measure Si {p) saturates at the 
value i. 



B. Measure based on quantum relative entropy 

Genoni et al [2l| have proposed, in a subsequent 
paper, a second measure of non-Gaussianity, this one 
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based on quantum relative entropy. They define non- 
Gaussianity of a state p as 

5^{p) = S{f)-S{p), (8.2) 

where S{-) is the von Neumann entropy of the state in 
question and f is the Gaussian state with the same first 
and second moments as the given state p. 

At first sight it would seem that 52{p) and our mea- 
sure A/'(/5) are very similar, the only difference being that 
p is replaced by Q{z) and that the trace operation in 
the formula for the von Neumann entropy is replaced in 
our measure by a phase space integral. A closer look 
reveals that this is not the case; (52 (p) does not reduce 
to M{p) under this kind of 'quantum-classical correspon- 
dence'. And 52{p) and N{p) turn out to be quite different 
entities. 

A qualitative difference between 52{p) and N{p) be- 
comes manifest when one compares these two measures 
in the context of a pure state. As the von Neumann 
entropy of a pure state is zero, (52 (/5) reduces to S{t), 
the von Neumann entropy of the Gaussian state with the 
same first and second moments as p. In other words (52 (p) 
does not consult, in the case of pure states, moments or 
cumulants of p of order higher than 2. Consequently, all 
pure states which have the same set of first and second 
moments but differ in higher moments will get assigned 
the same non-Gaussianity (52 (p). This is not the case with 
our measure N{p). 

To bring out a second qualitative difference we check 
if (52 (p) satisfies the shape criterion. To this end we ask 
if (52 (p) will ascribe the same amount of non-Gaussianity 
to the PATS and the photon number state, i.e.,, whether 
(52 (p) evaluated for the PATS Pp^^s independent of 
the temperature parameter x. We find that this is not 
the case. This is shown in Fig. ^ wherein we present 

(52(PpATs)' ^'^'^ fixed value m = 1, as a function of x. 

We conclude this Section with a further remark. With 
reference to Figs. ^ and (g]), while the non-Gaussianity 

measures ^i(Ppats) ^'^^ '^2(Pp™xs)' fixed m, vary with 
the temperature (or scale) parameter cc, thus failing the 
shape criterion, the variation is not monotone. The sig- 
nificance of the temperatures at which these measures 
assume their respective minimum values is not clear. 

IX. CONCLUDING REMARKS 

We have presented a measure of non-Gaussianity of 
quantum states based on the Q function. In doing 
so we have been guided by the fundamental principle 



that any measure of non-Gaussianity is an attempt to 
make a quantitative statement on the departure of the 
shape of the Q function from Gaussian, and the mea- 
sure must therefore remain invariant under all transfor- 
mations which do not change the shape of the Q function. 

Uniform scaling of all the phase space coordinates at 
the level of the Q function has proved to be an important 
shape preserving transformation, and our shape criterion 
demands that non-Gaussianity of the photon-added ther- 
mal states should be independent of temperature. 

We have explored various properties our measure 
which meets the shape criterion. We have presented ana- 
lytical and numerical results on the non-Gaussianity of a 
few families of quantum states. We have also compared 
our measure with other measures of non-Gaussianity 
available in the literature. 

Our measure M{p) meets the shape criterion which, in 
our opinion, should be respected by every good measure 
of non-Gaussianity. We hasten to add, however, that this 
is not the only measure that meets this criterion. For 
instance, if --yf^") is an appropriate linear combination 
of the cumulants of order 2n, and 7^^^ an appropriate 
linear combination of the cumulants of order 2, it is clear 
that the ratio between 7^^"^ and the n*"^ power of 7'^' 
will meet this criterion, for every n > 2. Our choice 
N{p) has the attraction of being immediately related to 
well-known entities like the Wehrl entropy and KuUback- 
Leibler distance. 

In the case of classical probability defined on a 2n di- 
mensional space C", one would have required the non- 
Gaussianity measure to be invariant under the full Eu- 
lidean group consisting of tanslations and all S0(2n) ro- 
tations. In the case of phase space, S0(2n) rotations 
which fall outside the subgroup Sp(2n, 7^) n S0(2n) are 
unphysical, and hence the restriction to this subgroup of 
passive linear systems. 

Our shape criterion rests on the invariance semi-group 
of Q functions which is different from the invariance semi- 
group of the Gaussian family of states - operations which 
map Gaussian states into Gaussians. The latter semi- 
group includes the full Sp(2n, TZ), and not just the inter- 
section subgroup Sp(2n, TV) n S0(2n). It further includes 
a whole family of completely possitive maps known as 
Gaussian channels. 

Finally, it may be noted that a scaling transformar- 
ion similar to the one we have implemented on Q func- 
tions cannot be implemented on Wigner functions. This 
follows from the following fact: that W{a) is a Wigner 
function docs not imply that A^" W{\a) is necessarily a 
Wigner function. 
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